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exact quantum dynamics and the quantum phase model 
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A simple four-mode Bose-Hubbard model with intrinsic time scale separation can be considered 
as a paradigm for mesoscopic quantum systems in thermal contact. In our previous work we showed 
that in addition to coherent particle exchange, a novel slow collective excitation can be identified by 
a series of Holstein-Primakoff transformations. This resonant energy exchange mode is not predicted 
by linear Bogoliubov theory, and its frequency is sensitive to interactions among Bogoliubov quasi- 
particles; it may be referred to as a second Josephson oscillation, in analogy to the second sound 
mode of liquid Helium II. In this paper we will explore this system beyond the Gross-Pitaevskii 
mean field regime. We directly compare the classical mean field dynamics to the exact full quantum 
many-particle dynamics and show good agreement over a large range of the system parameters. 
The second Josephson frequency becomes imaginary for stronger interactions, however, indicating 
dynamical instability of the symmetric state. By means of a generalized quantum phase model for 
the full four-mode system, we then show that, in this regime, high-energy Bogoliubov quasiparticles 
tend to accumulate in one pair of sites, while the actual particles preferentially occupy the opposite 
pair. We interpret this as a simple model for thermal pressure. 
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I. INTRODUCTION 

Thermodynamics remains one of the most complete 
and powerful physical theories in our possession. The 
connection to statistical physics established in the 19 th 
century was a very effective first attempt to understand 
the microscopical foundations of thermodynamics and to 
combine this highly effective theory with classical me- 
chanics. Yet even today the exact relationship between 
thermodynamics and microphysics is not completely un- 
derstood. Contemporary experiments investigating small 
quantum systems, for example Bosc-Einstein condensates 
in optical lattices, are a recent advance that revives old 
questions. The process of equilibration and thcrmaliza- 
tion in such small and isolated systems has received con- 
siderable theoretical re-examination [1-5] , and motivated 
the proposal of new concepts of equilibrium [6-11]. 

The approach we wish to pursue in this context is based 
on the view that thermodynamics is ultimately an adi- 
abatic theory, in the dynamical sense of deriving effec- 
tive Hamiltonians by eliminating high-frequency degrees 
of freedom. This adiabatic basis of the subject is often 
concealed, in statistical mechanics, by the definition of 
ensembles; but the fundamental justification of ensem- 
bles lies only in representing time-averaging. For small 
systems, however, it may become feasible to return to 
the Hamiltonian and apply adiabatic analysis directly. 
Such an approach may then shed light on the actual mi- 
crophysics underlying thermodynamics, and on its basic 
relationship with quantum mechanics. We therefore be- 
gin this ambitious program for quantum thermodynamics 



with its simplest possible example, and consider a mini- 
mal model Hamiltonian for two small quantum systems 
in thermal contact. 

At least if we restrict attention to models that some- 
what resemble the confined gases to which classical ther- 
modynamics traditionally refers, the smallest nontrivial 
quantum system is surely the two-mode Bosc-Hubbard 
(BH) model. This simple model has been extensively 
studied, not only theoretically [12-15], but also exper- 
imentally, since it is accessible in many cold atom lab- 
oratories [16-19]. Weakly coupling two of these simple 
systems together might therefore be the simplest case of 
two quantum systems in contact, such that they can ex- 
change both energy and particles independently. In our 
previous work [20] we introduced a four-mode BH system 
with a clear time scale separation as such an example of 
two exactly similar quantum subsystems in thermal con- 
tact, and showed how something similar to heat exchange 
may be described in familiar quantum mechanical terms. 

Each subsystem in this model is itself two-mode BH 
model - a pair of bosonic modes coupled so as to repre- 
sent an idealized Josephson junction (J J) [13, 14, 21-23]. 
We specify that the subsystems' characteristic timescale- 
their Josephson frequency-is to be the shortest one in 
the entire system, so that observation and control can 
address only much longer time scales. This leads to a 
Hamiltonian of the form 

H = H L + H R + H LR 

H hR = - | (V 1L a 1R + 4l«2R + H - c) , 
H a = (a\ a a 2a + a 2 > la ) + U^af a a ia , (1) 
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with on-site interaction U and index a = L, R for the left 
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and right subsystem, respectively. Coupling between the 
two subsystems is modeled by an additional Josephson 
tunneling between them with a much smaller rate coeffi- 
cient ijj than that of the intra-subsystem tunneling £1 to 
implement time scale separation (uj/tt <C 1). 

In [20] we showed that although only the total particle 
number N is an exactly conserved quantum number, in 
adiabatic theory there exists another conserved quantity, 
the total number of Bogoliubov exciatations J of both 
of the subsystems. We called these excitations 'josons', 
since they are the elementary excitations of Josephson 
junctions. We showed that although linear Bogoliubov 
theory identifies only one low-frequency normal mode in 
this system, in fact there exists a second low-frequency 
collective mode in this adiabatic limit, which emerges 
from the beating of the two high-frequency Bogoliubov 
modes, and corresponds to the exchange of josons be- 
tween the subsystems. Thus, the two subsystems can 
exchange particles by ordinary ('first') Josephson oscilla- 
tions on a 'slow' time scale, but they can also exchange 
josons (i.e. heat, in the sense of high-frequency energy), 
as the amplitudes of their fast Josephson oscillations beat 
slowly back and forth. We call this second low-frequency 
collective mode 'second Josephson oscillations', by anal- 
ogy with first and second sound in liquid Helium II. 

While the second Josephson mode is not an elementary 
excitation within the full theory, since it does not exist 
as a perturbation of the zero-joson ground state, it does 
appear as an elementary excitation within the adiabatic 
effective theory that describes low-frequency excitations 
around highly excited stationary states. Moreover, post- 
Bogoliubov nonlinear effects that are resummed in the 
adiabatic analysis shift the second Josephson frequency 
away from the simple Bogoliubov beat frequency, giv- 
ing it a nontrivial collective character. Within the low- 
frequency effective theory, however, the two Josephson 
modes are described in exactly the same way, with a low- 
frequency Hilbert space spanned by quantized excitations 
of both. Thus, our combined four-mode BH system may 
be regarded as a simple toy model for quantum heat ex- 
change. 

Our analysis of these second Josephson oscillations in- 
cluded the derivation of an extended free quasiparticle 
theory for collective excitations of the system consid- 
ered, simply by generalizing the Bogoliubov quasiparticle 
to treat adiabatically invariant quasiparticle numbers as 
fixed. Despite the fact that our derivations were made 
in terms of operators, the whole procedure was in princi- 
ple classical, since we restricted ourselves to the regime 
of large N and the limit J <C N. Thus our results in 
[20] could be presented using nonlinear Gross-Pitaevskii 
mean field simulations. 

In the present paper we now adopt two different ap- 
proaches to go beyond this mean field regime. In the 
first place we extend our numerical analysis to the exact 
full quantum dynamics of the four-mode system. The to- 
tal number of particles N is restricted by computational 
effort, but it is achievable to solve cases with TV > 40. 



While the large N theory of [20] requires corrections for 
small N or for J ~ N (presented elsewhere [24] ) , we show 
here that the quantum results do indeed agree with the 
mean field approximations for larger N. 

Secondly, then, we extend the standard quantum phase 
model (QPM) for Josephson junctions [25-29] to the two 
coupled junctions of the four-mode system. This can 
easily be achieved by requantizing the classical canonical 
conjugate variables. This model, being also only valid 
for large N (although it might be generalized for all val- 
ues of N in the spirit of [30]), nevertheless is able to 
predict small quantum corrections to the mean field so- 
lutions, such as quantum depletion. The QPM can also 
conveniently describe so-called 'running' states, in which 
josons preferentially are in one half of the system. We 
will show that in this regime the concentrated josons ef- 
fectively repel the atoms, pushing a proportion of them 
into the other half of the system. We interpret this effect 
thermodynamically in terms of josons providing pressure, 
similar to the thermal pressure of a macroscopic ideal gas. 



II. REVIEW OF THE JOSON THEORY FOR 
THE FOUR-MODE SYSTEM 



As a first step we implement a number-conserving Bo- 
goliubov theory for the two subsystems. To that end 
we define atom- moving operators a^ a ,a a by a Holstein- 
Primakoff transformation (HPT) [31] of the form 

2\J N a - ala a a a = (a\ a + a\ a ){a la - a 2a ). (2) 

These new operators commute with N a = J2i=i 2 ^la^ia 
and fulfill the usual bosonic commutation relations 
[a a , fit ] = 1. Thereupon we perform a Bogoliubov trans- 
formation by letting a a = u a b a + v a tf a and choosing 
u ai v a so as to diagonalize the quadratic terms of H a 
in (1) while demanding [b Q ,&|J = 1. Being interested 
in dynamics slow compared to the timcscale set by 0, 
i.e. an adiabatic regime, we now apply a rotating wave 
approximation (RWA) and ignore all terms that do not 
commute with the leading term proportional to w a b a in 
H a . In these terms the subsystems' Hamiltonian reads 



n 

2" a ' 2 
U m + 2UN, 



K + ^N a (N a 2) + y/n(n + 2UN a ) St S Q 
blX + OiUN- 1 ). (3) 



4 Q + 2UN a 

Note that the number-conserving Bogoliubov formalism 
introduced here via an HPT is somewhat different from 
the standard number-conserving Bogoliubov theory [32- 
35] , but follows quite naturally from the properties of the 
HPT. 

In the new variables the Hamiltonian of the subsystems 
(3) is diagonal. The coefficient of the quadratic term is 
the well known Josephson frequency of the excitations 
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created by b\, which we call 'josons', since the are the 
elementary excitations of the two Josephson junctions. 
The coefficient of the term quartic in b a is in general of 
the same order as U, but has opposite sign. This im- 
plies the simple but somewhat surprising fact that when 
atoms repel each other, josons attract each other; and 
vice versa. 

The total Hamiltonian of the whole four-mode system 
in zeroth order in oj/Q, is just equal to £1J, where Q — 
+ UN) is one of the high Bogoliubov frequencies 
of the four-mode systen and 



J 



(4) 



the total joson number. In the terms of Hamiltonian 
adiabatic theory J is an adiabatic invariant. Thus by 
maintaining our previous RWA we can now compute low 
frequency dynamics to leading order in small frequency 
ratios, dropping all terms in H that do not commute 
with J. The total atom numbers on each side, however, 
are not separately conserved for nonzero coupling oj. Ex- 
pressing the coupling term in terms of the Bogoliubov 
modes therefore forces us to again define new bosonic 
annihilation and creation operators A a , A* a that fulfill 
A' a A a = N a . These operators, which commute with the 
previously defined aj,, a a (and therefore naturally also 
with the Bogoliubov transformed quantities) enable us 
to take into account the rearrangement of atoms between 
the two subsystems. In these terms the coupling Hamil- 



tonian H hR assumes the form 



H, 



LR 




A T 



«R«R 

A D 



(5) 



+ 



^l^rA^r 




where by applying the same RWA as before we may write 
ala a = (ul + v 2 a p a b a , a[a R = UlUrS^Sr + v L v R b R b L . 

Assuming that our system is in an eigenstate of the con- 
served total atom number N = N L + N R with eigenvalue 
N, and of the total joson number J with eigenvalue J, 
we can now perform two final HPT. One to express N a in 
terms of TV and an A-conserving operator a that trans- 
fers atoms between the subsystems defined by 

A L R = 1 -[N± tfV V + a)} + 0{N-W) (6) 

and a completely analogous one to express b a in terms 
of J and a J-conserving operator b that transfers josons 
between the subsystems defined by 



Jl,r - \[J ± J 1/2 (& f + b)] + 0(J- 1/2 ). 



(7) 



To identify collective excitations, we finally linearize H 
in a and b to obtain for fixed values of N, J 

H lin = totfa + jN{a) + a) 2 + wjVb + ^ J(6 f + bf 
+ ~(Njy/ 2 (al + a)(tf+b), (8) 



with the definitions 

n + UN/2 



to j = oj 



^0(0 + UN) ' 



u., =- 



[/ 4f2 + UN 
2" fl + UN 



(9) 



The first line of (8) shows two different forms of the 
standard single JJ Hamiltonian: one providing Joseph- 
son oscillations of atoms between the L and R su bsys- 
tems, with Josephson frequency Q = y / w(o7+ UN); and 
another one implying Josephson oscillations of josons, 
i.e. seco nd Josephson o scillations, with the frequency 
ujj — y/ojj(ojj + UjJ). These two different types of 
Josephson modes are exactly analogous, except that, as 
mentioned before, Uj and U are of opposite sign. The 
second line of (8) finally couples the atom and joson 
modes, and arises because the frequency of fast Joseph- 
son oscillations in each subsystem depends on the respec- 
tive atom number. 



III. FULL QUANTUM DYNAMICS 

The crucial test of the joson theory will now be com- 
parison to actual quantum dynamics, in particular in re- 
gard to the frequencies of collective excitations. The two 
modes appearing in the final Hamiltonian (8) derived in 
the previous section are still coupled, but they may eas- 
ily be decoupled to yield the collective mode frequencies 
predicted by our analytical adiabatic theory: [20] 



lo 2 + uj 2 , 



± 



ujujjW 2 NJ 
n + UN 



"I 1/2 



(10) 



We will then test our adiabatic approximations by com- 
paring this prediction with the results of numerical exact 
diagonalization of the full four-mode Hamiltonian. Since 
the dimension of the Hilbert space is only Dim('H) = 
0(N 3 ), with destktop computers we may handle total 
particle numbers of more than 40. Because J is an adi- 
abatic invariant, though not an exactly sharp quantum 
number, the expectation values of J in exact eigenstates 
of the Hamiltonian are all very close to integers, and 
possess very small quantum uncertainties. Therefore we 
may unambiguously assign (N, J) labels to every eigen- 
state by setting J to be the rounded value of (J). To 
numerically check our joson theory we then take a super- 
position of the two lowest energy eigenstates for given 
values of N and J as initial state, and evolve it in time. 
Then as a function of time we plot the expectation value 
of the atom difference AA = (A L — A R ) and the joson 
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difference A J = — ^r^r)- A typical example of 

the resulting evolution is shown, or N — 44 particles, in 
Fig. 1. We observe essentially perfect agreement of our 
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FIG. 1. Time evolution of AJ = <S£& L - b R 6 R ) and A7V = 
{N L — iV R } of the full quantum system for the parameters 
A = 1, u = 0.1, UN = 0.1, N = 44 and J = 5. Red +: 
numerical A J. Red solid line: oscillation at second Josephson 
frequency w_ computed in the text, fitted for phase. Black 
o: oscillation at beat-frequency ojj. Blue X: numerical AN. 
Blue solid line: oscillation at a beat of the frequencies to + and 
w_ fitted for amplitude and phase. 

analytically derived collective mode frequencies with the 
numerical results. This is not trivial, since what we have 
derived amounts to a free quasiparticle theory for excita- 
tions around higher excited states of our system, rather 
than around the ground state. The josons indeed oscillate 
between the subsystems with second Josephson frequency 
cD_ , which is the lowest frequency in the system, and in 
particular is shifted well below the Bogoliubov beat fre- 
quency u>j (shown in our plots for comparison), proving 
that the second Josephson oscillation is sensitive to post- 
Bogoliubov nonlinear dynamics. Since atoms and josons 
are coupled, our initial state also includes an excitation 
of the first Josephson oscillation of the atoms. 

For another comparison we can also show the results 
of mean field evolution of our system, under a four-mode 
Gross-Pitaevskii nonlinear Schrodinger equation. In con- 
trast to the dynamical procedure used in [20] to excite a 
mean field state with some value of J, we now directly 
translate a given joson number to a certain particle im- 
balance, by inverting our series of Bogoliubov transfor- 
mations, in order to construct an initial mean field state 
with an exact value of J. The resulting mean field dy- 
namics for the same parameters as in the quantum case is 
shown in Fig. 2. Here again, the agreement with the lin- 
earized joson theory is excellent, as already shown in [20]. 
We may also confirm that for N — 44 particles the agree- 
ment between mean field theory and the exact many par- 
ticle dynamics is already very good. To explore instead 
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FIG. 2. Time evolution of AJ and AN of the corresponding 
mean field system for exactly the same parameters as used in 
Fig. 1. 



the less classical regime of small N, while keeping large 
enough J to have josons present to oscillate, the condi- 
tion J <C N must be abandoned. This requires relaxing 
several simplifying assumptions made in the derivation 
in section II, in order to construct a more general theory 
of collective excitations including josons, which will be 
valid for small N and for J ~ N [24]. 



IV. FOUR-MODE QUANTUM PHASE MODEL 

The standard procedure of deriving a QPM for a single 
JJ starts with the two-mode Gross-Pitaevskii equation 
for the classical (mean field) variables, which is obtained 
from the quantum mechanical Heisenberg equation of 
motion by replacing all the operators with c-numbers. 
The procedure is to represent this classical theory in po- 
lar canonical variables, obtaining a Hamiltonian func- 
tion with a cosine-potential. Re-quantizing canonically 
then leads to the Hamiltonian operator of the QPM, for 
which the corresponding Schrodinger equation is a Math- 
ieu equation. Because canonical quantization in polar co- 
ordinates is not generally valid, the QPM Hamiltonian 
operator is an approximation, and is only accurate in the 
large N limit. Unlike the mean field theory on which it 
was based, however, the QPM can express quantum su- 
perpositions, including the Gaussian squeezing that rep- 
resents quantum depletion, as well as highly non-classical 
'running states'. 

We therefore now construct a QPM for our four-mode 
system, which will respect the existence and conservation 
of josons as well as of atoms. We start with the Hamilto- 
nian obtained above after the first HPT and Bogoliubov 
transformation, given by the sum of (3) and (5). The cor- 
responding mean field Hamiltonian function is obtained 
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by replacing the operators according to 



(11) 



for a = L, R. In terms of these new canonically conjugate 
pairs (Ha, flu) and (r a ,f3 a ) the Hamiltonian function be- 
comes 

- kVr L r R u L u R cos(i? R - i? L + /3r - (12) 



+ v h v R cos($ R - t? L - /3r + /?l) 

- w [(i? L - K + «£)rL)(i2 R - (4 + v 2 R )r R )] 
x cos(tf R - 1? L ), 



1/2 



omitting constant terms. As in the standard procedure 
for a single JJ we now introduce relative populations 
R = (i?L — Rr)/2 = AN/2 and r = (r L - r R )/2 = A J/2 
as well as relative phases i9 = $l — $r and /3 = /?l — /? R 
of atoms and josons as new conjugate variables. Expand- 
ing to second order in R and r, and assuming the high 
occupancy limit for atoms and josons, the Hamiltonian 
function may be cast into the form 

"Hqpm = UR 2 + Ujr 2 + 2U^lrR -^{N-J) cos?? 

— ^ J cos i? cos/? + ^ Jsin^sin/3 (13) 

where we have used the abbreviations luj and Uj as 
introduced in (9), and defined the new quantity J = 
(CI + UN/2)J/Vl. We then follow the standard proce- 
dure to re-quantize "Hqpm canonically. Since R and /3, 
r are pairs of canonically conjugate variables, the QPM 
approach is to do this simply by replacing them with 
operators expressed as 

# = R = -id#, j3 = P, r = -\d p . (14) 

This leads to the Schrodinger equation id t ^ = Hq PM ^ 
for the wave function ^(d,/3). Since d and (3 are phases, 
^(i?, /3) must be 27r-periodic in both arguments. In the 
framework of this QPM we are now able to analyze not 
only Josephson states, with the support of localized in 
the potential wells around $ = (3 = 0, but also highly 
excited running states. 

To calculate the frequencies of excitations around the 
most localized Josephson states, we may simply expand 
the cosine-potentials up to second order and diagonalize 
the resulting coupled oscillator Hamiltonian. One of the 
two resulting frequencies a/X±, however, is imaginary, 
since A_ < while still A + > 0. This may be under- 
stood easily, since explicitly performing the diagonaliza- 
tion shows that these frequencies indeed coincide with the 
quantum phase model limit J, N — >• oo of the frequencies 
(10), which can be written as Q± = uj 2 + \± up to higher 
order corrections. We already know from the linearized 



theory reviewed above that because the inter-joson and 
inter-atom interactions are opposite in sign, either the 
first or second Josephson mode will become dynamically 
unstable for high enough J or N. What linearization 
of the QPM around the phase-localized Josephson states 
has now shown is that system must already have reached 
this dynamical instability threshold when the QPM limit 
applies. The standard QPM can therefore never describe 
both forms of Josephson oscillation. 

We may nevertheless usefully apply the QPM to our 
system, to describe the non-phase-localized states that do 
become stable in the large J, N regime. For simplicity we 
consider only the case of repulsive inter-atomic interac- 
tions, with attraction between josons. With Uj < 0, the 
energy decreases with f 2 , so low energy favors the largest 
possible eigenvalues f 2 . These eigenvalues can be as large 
as J 2 , and so in the large J limit where the QPM becomes 
accurate, the Ujf 2 term in (13) will always dominate the 
low-energy states. We therefore treat the coupling terms 
between josons and atoms in (13) as a perturbation, and 
assume that to zcroth order the wave function ^ may 
be decomposed as the product of a running state for the 
josons and a phase-localized state for the atoms, of the 
form 



1 



2ir 



(15) 



Here r is the eigenvalue of f, and the low energy states for 
fixed J will have |r| ~ J >• 1. These J 'running' eigen- 
states of f represent states with many high-frequency jo- 
son excitations, in which the majority of josons remain 
permanently localized in either the L or R half of the 
four- mode system, depending on the sign of r. This cor- 
responds to insulator-like quantum number locking in the 
second Josephson mode, while the first Josephson mode 
retains its usual phase-localized behavior. 

In a Born-Oppenheimer sense we may consider r to be 
sharp, and apply an r-dependent adiabatic Hamiltonian 
H r for the first Josephson mode alone: 



/ 

/ rr 
JO 



(fl,f3)H QPM y rn ($,f3)df3 



(16) 



Ug 2 - ^{N - J) costf + r 2 ( Uj - U^- 

2 \ n 2 



where we have defined the shifted atom imbalance g = 
—ids + rfl/fl. Thus the effective Schrodinger equation 
for the first Josephson mode, in a number-locked state of 
josons with given r, is again the Mathieu equation of the 
standard QPM, with the eigenstates given by Mathieu 
functions. There is a nontrivial effect from the josons on 
the atoms, however, in that the atom population imbal- 
ance AN = 2(R) is shifted away from zero by an amount 
proportional to the the joson imbalance r: 

AN = 2(g)-2r2=2(g)-Aj2. (17) 
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This implies that if the mutually attractive josons are 
concentrated in the L subsystem, then the mutually re- 
pulsive atoms are partially driven toward the R subsys- 
tem; and vice versa. 

From this we conclude that josons exert something like 
pressure on the atoms. This further supports our tenta- 
tive, toy-model interpretation of josons as quanta of heat, 
since a hot gas should exhibit thermal pressure. Col- 
lecting josons on one side means heating this subsystem, 
according to our dynamical toy model for microscopic 
thermodynamics as an adiabatic theory. It seems con- 
sistent with our picture, therefore, that the hotter gas 
should reveal its higher pressure, by forcing more atoms 
into the other subsystem. 

Of course, given that a certain number of high- 
frequency joson quanta are localized in one half of the 
system, it is obvious that the system can lower its energy 
by having fewer atoms present in the half of the system 
where they will have to participate in high-frequency ex- 
citations. Precisely because this is such a simple effect, 
it is not clear that it is entirely an artifact of our toy 
model. It may perhaps be a fair representation of the mi- 
crophysics of thermal pressure in more general systems. 

V. DISCUSSION 

We have argued that time scale separation is the essen- 
tial dynamical element in hydrodynamic collective exci- 
tations of many-body systems. Such excitations are very 
common; they include ordinary sound waves in air. In 
macroscopic cases they may be described accurately by 
a low-frequency effective theory in terms of thermody- 
namical fields like pressure and temperature. The most 
natural relationship to postulate between these classical 
thermodynamical quantities, and the quantum Hamilto- 
nian microphysics that underlies them, is that the ther- 
modynamic quantities are approximate adiabatic invari- 
ants in the full system dynamics. They are not changed 
by rapid microscopic motion, but only by much slower 
macroscopic kinetics, via processes like convection or dif- 
fusion. The way in which the simple description of these 
adiabatic degrees of freedom emerges is certainly not sim- 
ple, however. Even when the low-frequency effective hy- 
drodynamics may be accurately linearized, as for small- 
amplitude sound waves, this linear theory will typically 
represent some strong and complicated nonlinear effects 
at the microscopic level, for example in the collisions be- 
tween particles that maintain local equilibrium. 

With our four-mode BH model we have explicitly real- 
ized this general scenario for relating thermo- and hydro- 
dynamics to Hamiltonian microphysics, by constructing 
a simple quantum many-body Hamiltonian with a clear 
time scale hierarchy Simple as it is, this model nonethe- 



less has a quite complicated spectrum and exhibits quite 
nontrivial dynamics. We were able to show that there 
exists a highly excited collective mode that represents 
the exchange of high-frequency energy (josons) between 
the two subsystems, with the exchange proceeding as an 
oscillation at low frequency. 

In this paper we extended the analysis of this nontriv- 
ial collective mode beyond our previous mean field limit, 
and we have done so in two ways. First we extended 
our numerical studies to the quantum case, and showed 
excellent agreement of our linearized joson theory with 
the actual quantum evolution, for moderately large num- 
bers. We thereby confirmed that our second Josephson 
frequency directly manifests itself in the quantum spec- 
trum of the complete four-mode system, as an excitation 
energy above the lowest energy eigenstate with fixed iV 
and J. To extend our adiabatic analysis further into the 
quantum regime than this, it is necessary to relax the 
condition of large N and J <C N that we assumed intro- 
ducing the joson theory in [20] and also used here. The 
extended theory, valid for small N and J N, is described 
elsewhere [24]. 

Secondly we applied the quantum phase model to our 
four-mode system, and by examining number-locked jo- 
son states represented by running states of the quantum 
phase model we were able to show that localized josons 
exert a kind of pressure on atoms, in keeping with their 
postulated role as quanta of heat. Of course, the fact 
that josons concentrate spontaneously, in the case of re- 
pulsive interactions among atoms, does not represent the 
usual diffusive behavior of heat. The attractive inter- 
actions of josons need not entirely rule out their heat 
interpretation, however, since it is well known that spe- 
cial kinds of dynamics can produce systems with nega- 
tive specific heat capacity, e.g. self-gravitating systems, 
for which heat does indeed spontaneously concentrate, 
without any violation of thermodynamic laws. 

For further tests of dynamical conjectures about ther- 
modynamics, it will clearly be necessary to investigate 
larger systems with more degrees of freedom. In larger 
Bose-Hubbard systems there are, on the one hand, many 
new possibilities for the dynamics of josons. On the 
other hand, however, they in general tend to be dynam- 
ically chaotic, which might play an important role in 
the emergence of thermodynamics. If at some point the 
complexity of these models becomes theoretically com- 
pletely intractable, it may still prove possible to inves- 
tigate them by means of experiments on cold, trapped 
Bose gases. While such experiments would be difficult, 
they are clearly feasible in principle as an extension of 
demonstrated capabilities. They would be a fine exam- 
ple of the principle that we do not study physics to learn 
about cold atoms, but rather study cold atoms to learn 
about physics. 
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